Long time dynamics following a quench in an integrable quantum spin chain: 
local versus non-local operators and effective thermal behavior 
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We study the dynamics of the quantum Ising chain following a zero-temperature quench of the 
transverse field strength. Focusing on the behavior of two-point spin correlation functions, we show 
that the correlators of the order parameter display an effective asymptotic thermal behavior, i.e., 
they decay exponentially to zero, with a phase coherence rate and a correlation length dictated by 
the equilibrium law with an effective temperature set by the energy of the initial state. On the 
contrary, the two-point correlation functions of the transverse magnetization or the density-of-kinks 
operator decay as a power-law and do not exhibit thermal behavior. We argue that the different 
behavior is linked to the locality of the corresponding operator with respect to the quasi-particles of 
the model: non-local operators, such as the order parameter, behave thermally, while local ones do 
not. We study which features of the two-point correlators are a consequence of the integrability of 
the model by analizing their robustness with respect to a sufficiently strong integrability-breaking 
term. 

PACS numbers: 75.40.Gb, 75.10.Pq, 73.43.Nq, 03.65.Sq 



I. INTRODUCTION 

The study of the non-equilibrium dynamics of strongly 
correlated quantum many body systems has lately re- 
ceived an increasing amount of attention. While the 
first theoretical studies in this area were performed a few 
decades ago 1 - - — , interest on this subject has been confined 
for a long time to the theoretical literature. The main 
trigger behind the recent experimental advances on this 
subject has been the impressive progress in manipulat- 
ing cold atomic gases, which resulted in the first exper- 
iments exploring the coherent non-equilibrium dynamics 
of strongly interacting quantum systems in a controllable 
way£~-. 

On the practical side, the main advantage of studying 
many-body physics with cold atoms is the detailed mi- 
croscopic knowledge of the Hamiltonian describing these 
systems. This fact, together with the possibility of con- 
trolling the system's parameters with high accuracy and 
with the availability of long coherence times, has made it 
possible to study the unitary dynamics of many-body sys- 
tems 9 -^, such as Hubbard-like Hamiltonians 11 or artifi- 
cial spin chain model o 12 i 13 , a topic that can hardly be ad- 
dressed in traditional solid-state environments. Coherent 
non-equilibrium dynamics has been brilliantly demon- 
strated by the observation of collapse and revival cycles 
in systems quenched across a supcrfluid-to-Mott insula- 
tor quantum phase transition^ and in the study of the for- 
mation of topological defects during a quench of trapped 
atomic gases through a critical point££ 

One of the most fundamental topics in quantum sta- 
tistical mechanics which is presently being studied in 
this context is the connection between ergodicity, non- 



integrability, and thermalization in the dynamics of 
strongly interacting (but isolated) many-body systems. 
An important experiment addressing this issue was re- 
cently reported in Ref. [8j, where the absence of ther- 
malization in the coherent non-equilibrium evolution of 
a closely integrable one-dimensional Bose gas was ob- 
served. Motivated by these experimental findings, an in- 
creasing number of theoretical studies focusing on the 
dynamics following a sudden perturbation (a quantum 
quench) have been performed 1 ^—. 

Integrability is believed to play a crucial role in the 
relaxation process: in analogy to the well-known Fermi- 
Pasta-Ulam scenario in classical systems, integrable sys- 
tems are not expected to thermalize, but to be sensitive 
to the specifics of the initial stated— . This understand- 
ing was distilled into the proposal of a generalized Gibbs 
ensemble, keeping track of the initial value of all the con- 
stants of motion^ 8 - and constructed to describe the steady 
state reached after a quench. Several works have tested 
the conditions of applicability of such Gibbs distribution 
and its drawbacks 1 ^ -18 ' 21 ! 28-30 ! 34 . In particular, for a spe- 
cial quench in a ID Bosc-Hubbard model 2e », for integrable 
systems with free quasiparticles2£ and for the computa- 
tion of one-point correlation functions for a specific class 
of quench processes in otherwise generic integrable sys- 
tems^ 1 -, the long-time limit of the dynamics was shown 
to be well described by the generalized Gibbs ensemble. 
Moreover, the generalized Gibbs ensemble was shown to 
correctly predict the asymptotic momentum distribution 
functions for a variety of models and quenches 1 ^ 18 ' 29 . 
However, it should be pointed out that neglection of 
correlations of the occupation of different quasi-particle 
modes generally leads to incorrect predictions for the 
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noise and higher order correlators^. Turning to non- 
intcgrable systems, in this case thermalization is ex- 
pected to occur in general: this has been numerically ob- 
served in some circumstance s 19 ' 20 ' 22- — , while the transi- 
tion from integrable to non-integrable has been shown, at 
least in small lattices of interacting bosonic or fermionic 
particles 2 -, to take the form of a crossover. The mecha- 
nism of thermalization has been conjectured and numer- 
ically tested for certain system s 19 ' 20 to be analogous to 
the one proposed by Deutsche and Srednicki2£ for sys- 
tems with a classically chaotic counterpart. Nonetheless, 
it is worth pointing out that in some specific cases, like 
for gapped systems, the validity of this scenario has been 
questioned (the points raised include finite-size effects 2 ^ 
and the importance of rare events^ 9 -, which may drive the 
behavior of the long-time dynamics) , and the problem is 
still under debate. 

The purpose of this work is to go one step beyond 
the scenario proposed above and show that in quantum 
many-body systems the presence or absence of thermal 
behavior after a quantum quench does not exclusively de- 
pend on the integrability of the model, but also on the 
considered observable. In particular, we argue that, even 
in a completely integrable system, the asymptotics of the 
two-point correlation functions of an observable which is 
non-local with respect to the quasi-particle fields displays 
thermal behavior, while this is definitely not the case 
for a local one. This scenario may be violated in some 
specific and isolated cases, such as for the transverse- 
field correlations or the density-of-kinks operator after 
a quench towards the critical point. In this context, 
local and non-local operators refer to the structure of 
their matrix elements on the basis of quasi-particles: lo- 
cal means that the operator couples a finite number of 
states, while non-local means that it couples all states. 
To show the different behavior of the correlation func- 
tions of local and non-local operators, we study quantum 
quenches in a quantum Ising chain, focusing our atten- 
tion on the asymptotics of the two-point correlators of the 
order parameter o~ x and of the transverse magnetization 
a z after a quench of the transverse field. We explicitly 
show that the correlator of the order parameter, which 
is non-locally expressed in terms of the fermions which 
diagonalize the model, asymptotically displays a thermal- 
like behavior characterized by an exponential decay both 
in time and in space. As we already reported in Ref. [3lj |. 
the autocorrelation function of the order parameter after 
a quench takes the form 

WofWofiomve-W, (1) 

with a phase coherence time Tq depending only on the 
effective temperature T e g set by the ground state energy 
of the initial Hamiltonian \B) (below also called boundary 
state) , and on the energy gap A of the final Hamiltonian. 
Here we are going to show that such thermal behavior 
is not limited to the equal-site autocorrelation function, 
but also pertains to the large distance behavior of the 
equal-time correlator (asymptotically, for large times t), 



which behaves as 

(B\af +r (t)a*(t)\B)*e- r ^, (2) 

with a correlation length £q again determined by T e ff and 
A. Remarkably, as shown in Fig.QJa), the dependence of 
Tq on T ff (black circles) quantitatively agrees with that 
of the corresponding phase coherence time r^°- at ther- 
mal equilibrium as a function of the temperature (con- 
tinuous red line). A similar behavior is observed for the 
correlation length [see Fig. [2b)]. All of these results can 
be understood in terms of the generalization to the non- 
equilibrium case of a semiclassical analysis originally de- 
veloped for thermal equilibrium^ ., which explicitly takes 
into account the effect of quasiparticles generated after 
the quench [see, e.g., the dotted blue line in Fig.QJa)-(b)]. 
This thermal behavior is not present at all in the asymp- 
totics of the correlator of the transverse magnetization 
or in the nearest-neighbor spatial correlations of the lon- 
gitudinal magnetization (i.e., representing the density of 
kinks), that are operators which turn out to be local in 
the fermionic quasiparticles of the Ising chain. 

The fact that quasiparticles play a fundamental role 
in the physics of thermalization is not at all surpris- 
ing: the essence of quantum integrability is the possi- 
bility to describe a system in terms of quasiparticles sub- 
ject to purely elastic and factorizable scattering. A non- 
integrable system, in this sense, is just characterized by 
the absence of well defined quasiparticles: the scattering 
of excitations gcnerically leads to a cascade redistribut- 
ing the initial energy among the low-energy degrees of 
freedom. The results discussed in this paper indicate 
that if an observable is sensitive to the detailed scatter- 
ing properties of such degrees of freedom, i.e., it is local, 
then its asymptotics will be affected by the presence or 
absence of integrability. In turn, if such relation is non- 
local, then, as for the equilibrium physics, we expect that 
for its multi-point correlation functions only some gross 
low-energy features of the model are relevant, and de- 
termined by the universality class. We stress again that 
here locality is defined with respect to the quasiparticles, 
while the original degrees of freedom do not play a major 
role. 

The paper is organized as follows: in Sec. |TT]we define 
our model and the techniques used to study its dynamics 
after a quench. In order to clarify the meaning of locality 
in this context, we also discuss the continuum formula- 
tion of the model in terms of free Majorana fermions, 
and enlighten the physical content of the initial state \B) 
and its connection with boundary states in statistical 
field theory^. In Sec. Mil we summarize the main for- 
mulas for studying the time evolution of the system. In 
Sec. II VI we discuss the order parameter correlation func- 
tions: we consider in details the behavior after a quench 
showing that it is asymptotically thermal, both in time 
and in space. In Sec. [V] we focus on two examples of 
local operators (namely, the transverse field correlations 
and the density of kinks, which is essentially the nearest- 
neighbor space correlator of the order parameter) , which, 
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In the following we will consider sudden quenches of 
the transverse magnetic field strength T: after initial- 
izing the system in the ground state \4>(T )) = \B) of 
the Hamiltonian T-L(Tq), the field strength is suddenly 
changed at time t — from Tq to a new value T =/= Tq. 
Consequently, the state will evolve according to the new 
Hamiltonian 'H(r): 



m(r)t 



\B) 



(4) 
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Figure 1: (color online). Phase coherence time (a) and cor- 
relation length (b) extracted from the asymptotic decay of 
the two-point order parameter correlations, as a function of 
the effective temperature. Data are for quenches ending in the 
ferromagnetic phase of the Ising chain at V — 0.5 [see Eq. ©]. 
Symbols refer to the case of a quench and correspond to dif- 
ferent values of To (empty symbols are for To < T, while filled 
ones are for To > T); continuous red curves denote the equilib- 
rium values at finite temperatures, while dotted blue curves 
correspond to values obtained with a semiclassical analysis 
generalized to non-equilibrium cases. 



on the contrary, do not exhibit thermalization. The non- 
thermal behavior can be nonetheless described by means 
of a generalized Gibbs ensemble, as discussed in Sec. IV Cl 
A brief discussion of the effects of intcgrability break- 
ing is sketched in Sec. I VII where we provide numerical 
data showing the gradual disappearance of a typical non- 
thermal feature in the asymptotic spatial behavior of the 
order parameter correlators. Finally, in Sec. lVIIl we draw 
our conclusions. 



II. THE MODEL: GENERALITIES, 
CONTINUUM LIMIT AND BOUNDARY STATE 

In this paper we study the spin- 1/2 quantum Ising 
chain in a transverse magnetic field 4 ^, which is charac- 
terized by the Hamiltonian 



H(r) = -jV[ 



i=i 



Taj 



(3) 



where L is the number of spins in the chain and u" 
(a = x,y, z) are the Pauli matrices relative to the jth 
spin. Hereafter we impose periodic boundary conditions, 
and use a system of units J = 1, h = 1, and ks = 1- The 
quantum Ising chain is the prototype of an exactly solv- 
able quantum system and is characterized by two mutu- 
ally dual gapped phases, a quantum paramagnet (r > 1) 
and a ferromagnet (r < 1), separated by a quantum crit- 
ical point at T c = 1. 



Before entering into the details of the non-equilibrium 
dynamics of this model, let us briefly set the notation by 
reviewing the diagonalizatio n 43 i 44 of the Hamiltonian ([3]) . 
Introducing Jordan- Wigner fermions c\ defined by 



3-1 

Cj exp ( ITT 

1=1 



c\ci 



(5) 



and omitting constant terms, the Hamiltonian in Eq. ([3]) 
takes the quadratic form 



L-l 



u 



[4 c j+i + 44+1 + h.c. 



3 = 1 

+(-l) Np 



2r 7%i' / 

3 = 1 



c[ci 



44 



(6) 



The last term originates from the periodic boundary con- 
ditions imposed to the spins and its sign depends on the 
parity of the total number Np of c-fermions. Specifi- 
cally, if Np is odd, then all the bonds are identical and 
periodic boundary conditions on the fermions are im- 
posed (cl+i = ci). Antiperiodic boundary conditions 
(cl+i = — ci) are instead appropriate when Nf is even. 
The Hamiltonian in Eq. ^ conserves the fermion par- 
ity, therefore it can be formally split in two parts acting 
on different Fermion-parity subspaces, even (+) and odd 
(-): V. =H+ +n~, where = V ± UV ± denote the 
even/odd subspaces, and V ± the associated projectors. 
Since in the following we will consider initial states with 
an even number of fermions, we will focus our attention 
on the even sector only. 

The diagonalization of the Hamiltonian now proceeds 
in momentum space. Writing the c-fermions as Cj — 
Dfc e%hjc k, where k is 

fc^± 7r(2n r +1) with n = Q,--- ,L/2-l, (7) 

the Fourier representation of Eq. ([5]) in the even sector, 
H + , becomes a sum of independent terms 



(8) 



fe>0 
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where we introduced the Nambu vector cl = (c k , c_fe), 
and 

H fc = ( a h k - lbk ) (9) 
\ib k -a k J w 

with a k = 2(r — cos k) and b k — 2sinfc. This Hamilto- 
nian is easily diagonalized through a Bogoliubov rotation, 
defining the new fermionic variables 

A k = Rfc(r) c k (10) 

where A^ — (A k , A_ k )- The rotation R is explicitly 
given by 

R *< r > =(-«)■ 

Here 

£fe + a k ib k . 

u k = , = , v k = = , (12) 

y/2e k (e k + a k ) y/2e k {e k + a k ) 

while e k = sja\ + b\ = 2%/r 2 - 2rcosA: + 1 is the dis- 
persion of the quasiparticles in terms of which 

n+ = £ k(A{A k + Al k A- k - 1) . (13) 

fc>0 

The ground state of the system is the vacuum of the 
Bogoliubov quasiparticles defined in Eq. (fT0|). These be- 
come gapless at the quantum critical point r c , where the 
gap A = e = 2|1 — T| vanishes. 

B. Continuum formalism 

The continuum formulation of the quantum Ising 
model enlightens other properties, such as the local/non- 
local nature of its operators, and is helpful in understand- 
ing the properties of the initial state \B) associated to the 
quench process. Let us briefly discuss the main concepts 
of this formulation. 

In the vicinity of the quantum critical point, taking the 
continuum limit (a — > 0, where a is the lattice spacing), 
it is well known that the Ising model becomes a field 
theory of free relativistic Majorana fermions tp+ (r, t) and 
tp_(r,t)^. These are real fields, ip±(r,t) — ip±(r,t), and 
obey the equal-time anti-commutation relations 

{^(r,t),^(r',t)} = 6 tJ 6(r-r') . (14) 

The fermion particle excitations satisfy a relativistic dis- 
persion E(p) = a/ j? 2 + A 2 with a mass m = A, and their 
dynamics is described by the Hamiltonian 

% = J dr [itl) + d r ip+ — iip_d r ip- — 2irmp + ip^] . (15) 

The mode expansion of the fermion fields is given by 
r°° At- 

V>± (r, t)= — [a ± (k)A(k)e- iEt+ikr + h.c] (16) 

J — oo ^ 7r 



where {A(k),A^{k')} = 2%6{k-k'), and 

with uj — e l7r / 4 . In order to express later in a more concise 
form both the operators and the matrix elements of the 
model, it is convenient to re-write the mode expansion 
in a slightly different way. This can be done using the 
identity 




9{u) being the Heaviside theta, so that Eq. (|16p becomes 

iMr.t) - \/f / +0 ° ^e~ iEt +^ (S + i>(u) , 

V 2 J_ QO 2ttu 

(18) 

where E = A(u + u _1 )/2, p = A(u - u^ v )/2 and 

^{u) = ±=[6(u)A{u)-6{-u)A\-u)] . (19) 

The annihilation and creation operators, A(u) and 
A'(u), are the continuous version, in the variable u, of 
the lattice operators Ak and A\. The vacuum state |0) 
is identified by the conditions A{u) |0) = 0. Notice that 
ip(u) is an annihilation operator for u > and a creation 
operator for u < 0. It is also useful to define 

Mr,t) = \/f / + °° ^i,(u)e- iEt +^ . (20) 
V 2 J_ OQ 2ttu 

Since the work of Fradkin and Susskind^, it became well 
known that, in the continuum limit, the lattice operator 
of becomes 

[ ° ;f < : • (21) 

where a(r, t) is the "order parameter" , while fi(r, t) is the 
"disorder parameter" . On the other hand, for the lattice 
operator of we have 

/ e(r,t) foi'T > 1 
ar ^\-e(r,t) forr<l ' [ZZ> 

where e(r,t) is the "energy operator". Eqs. ([21]) and (|22|) 
express the self-duality of the Ising model. In particu- 
lar, for r > 1, <r{r,t) is a Zi odd operator with non-zero 
fermion number while /i(r, t) is a Z2 even operator with 
zero fermion number: so, a{r, t) has non-zero matrix el- 
ements only on odd number of fermions, while ^(r, t) 
has non-zero matrix elements only on even number of 
fermions. For T < 1, the role of the two operators is 
swapped and the situation is reversed: this symmetry is 
due to the self-duality of the model. These conclusions 
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are confirmed by the explicit expressions of the opera- 
tors in terms of the fermion fields given, for T > 1, by 
the normal ordered expressions^ 



e(r,t) = i : ip + (r, t) ip_(r, t) : 

/ i(r,t)=:e* t ) :, 

where p(r, t) is a quadratic form of fermions 

f du dv! (u — u') 
P{r,t) = -i \ — - 



(23) 



2nu 2txu' u + vf + iS 



(24) 



Notice that p(r, t) is expressed by a non-local quantity 
of the fermion modes which plays the role of the fermion 
string that accompanies the definition of the lattice op- 
erators (see Eq. (O). 

If we now set u — where (3 is the rapidity, using 
the operatorial expressions above, one can easily check 
that, for r > 1, the matrix elements of the various op- 
erators (between the vacuum state |0) and the asymp- 
totic states of n fermions identified by their rapidities 
|j8i, &,...,&,» are given by^8 



|e(0,0)|/?i,...,/3„) = 



iA sinh 




2n+l 



for n = 2 
otherwise 



a(0,0)\h,...,/3 2n+1 ) = (t) 2n+1 f[ tanh 



ft - ft 



i<j 



In 

(0| M (0, 0) 1)81, ... , ft>„) = « 2n f] tanh «^ . 

i<j 

These expressions of the matrix elements clearly show 
that the order/disorder operators are non-local fields 
with respect to the fermion excitations: they couple to 
an arbitrary odd/even number of quasiparticles and they 
may be consequently regarded as strongly-interacting 
fields, even though the Hamiltonian (fT5j) is quadratic. 
On the other hand, the energy field associated to a z is a 
local operator with respect to the fermions and it couples 
only to two-particle states. 

Another local field with respect to the fermion quasi- 
particle exitations is the density-of-kinks operator, whose 
lattice definition is 



(25) 



This operator is in fact at most quadratic in the 
fermion fields: this can be easily seen firstly us- 
ing the Wigner- Jordan fermions to express fjCj+i = 



c]c 



J+ 1 + c j c j+i 



h.c. 



and secondly rewriting this ex- 
pression in terms of the fermionic variables Ak , by means 
of the Bogoliubov transformation. Therefore this opera- 
tor only couples to two-particle states. 



C. Boundary state 

Before proceeding with the description of the dynam- 
ics following a quantum quench, it is worth discussing 
the representation of the initial state \B) in terms of the 
quasiparticles diagonalizing the final Hamiltonian and 
its connection to boundary states in statistical field the- 
ory^. Despite the fact that this connection can be eluci- 
dated directly within the lattice model^, here we follow 
an elegant route employing the continuum formulation 
and the equations of motion for the Majorana fields: 



(d t +d r )iP+(r,t) = 
(d t -d r )i>-(r,t) = 



(26) 



In this formulation the quench process consists of an 
abrupt change of the mass, Ao — > A, of the fermion 
field at t = 0. Since the equations of motion (|2"6")l are 
of the first order in <9 t , the fermion field should have no 
discontinuity at t = and therefore should satisfy the 
boundary condition 



V£M = o) = v±M = o) 



(27) 



where ip^. (r, t) and i[>± ( r i t) denote the fields relative to 
the masses A and A, respectively. Eq. (|27|) implies a 
linear relation between the modes of the field before and 
after t — 0. In order to compare directly with the lattice 
results, it is convenient to discretize the space variable in 
units of the lattice spacing a. In the following, focusing 
for simplicity on the case in which both To , T > 1 , let us 
write the mode expansion (|18[) of the two fermionic fields 
as 



(r,t)= / ^[a ± (k)A(k)e 
Jbz 2tt 



— iEt+ikr 



(28) 



where the integral is over the first Brillouin zone |fe| < 
w/a, E = A 2 +p 2 with p = (2/o) sin (ka/2), and 

a±{k) = w* 1 y Analogous expressions hold for the 



two components ip± (r, t) expressed in terms of the modes 
Ao(k) and Al(k), with E = ^A 2 +p 2 . 

Let us denote by \B) the ground state of the fermion 
fields ip±(r,t) and by |0) the ground state of the field 
"0±( r i^) : these states are annihilated by Ao(k) and 
A(k) respectively. Extracting at t = the Fourier 
mode ip±(k) of the fermion fields, defined by ip±{r,0) — 
f BZ dk/(27c)i, ± (k)e ikr , 



^+{k) = a + (k)A(k) + a+{-k)AH-k) 
■4>-(k) = a-(k)A(k) + aS{-k)A\-k) 



(29) 



(with analogous expressions for ip±) and imposing the 
boundary condition (|27p , one obtains the linear relations 
among the modes: 



A(k) = U(k)A (k) 
A\k) = U(k)A\{k) 



iV(k)Al(- 
iV(k)Ao(- 



k) 
k) 



(30) 
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Figure 2: With respect to the modes of the Hamiltonian at 
t > 0, the boundary state \B) appears as a coherent super- 
position of an infinite number of pairs of particles with equal 
and opposite momentum. 



where 
U{k) -- 
V(k) -- 



1 



1 



y/(E + p){E+p) + ^/{E -p)(E- p) 



2^EEo 



y/{E -p)(E+p) - y/(E +p)(E-p) 



Notice that U(k) = U(-k), V(k) = -V(-fc). The inverse 
relations of Eqs. (l30l) are given by 



(31) 



Ao(k) = U{k)A{k) + iV{k)A\-k) , 
Al(k) = U{k)AHk)-iV(k)A(-k) . 



Looking at these equations, we see that a quench in the 
mass of the fermion is simply equivalent to a Bogoliubov 
transformation of its modes - transformation that is ruled 
by the functions U{k) and V(k). The role of the initial 
state of the quench process is played by the ground state 
\B) of the fermion fields ijj±(r,t). 

It is simple to work out the properties of the state \B), 
hereafter called, for obvious reasons, boundary stat e - 17 > 41 . 
Since this is the vacuum of the fermion ijj°(r,t), we can 
use the formula 



(B\A (k)Al(k')\B) = 2w5(k-k') 



(32) 



and the Bogoliubov transformations (|30l) to easily com- 
pute the expectation values of modes A(k) and A'(k) on 
this state: 



{B\A(k')A(k)\B} 
{B\A{k')A\k)\B) 
(B\A\k')A(k)\B) 
{B\A\k')A\k)\B) 



2iriU(k')V{k')5{k' + k) 
2nU 2 {k')5(k' - k) 
2irV 2 (k')S{k' - k) 
2mU(k')V{k')6(k' + k) 



(33) 



However, we can gain new information on the nature of 
the boundary state \B) if we directly express it in terms 
the oscillators A(k) and A'(k). This can be done by 
using once again the condition that \B) is annihilated by 
all the operators Ao{p), but employing, this time, the 
first equation in (|31|) for the operators Ao(p). In this 
way, we get the following infinite set of equations for the 
boundary state \B): 

[U(k)A{k) + iV{k)A^{-k)] \B) = 0, (34) 



whose formal solution is a BCS-type state 

\B) = \{ [U{k) + iV(k)A\k)A\-k)\ |0) , (35) 

k 

where |0) is the vacuum of the fermions V-fc( r j0- Notice 
that, disregarding an overall normalization, the state \B) 
can also be expressed as 



\B) = exp 



BZ 



dk 
2^ 



JC{k)A\k)A\-k) 



|0> , (36) 



where fC(k) = U{k)/V(k) = -/C(-fc). Therefore this be- 
longs to the special class of initial states, recently studied 
in Ref. [2l|, represented by a coherent superposition of 
particle pairs for which the long-time limit of one-point 
functions is described by a generalized Gibbs ensemble. 
The expression (|35l) or, equivalent ly ([36]) . shows that, as 
a state of the Hilbert space of the fermion ip(r,t), the 
boundary state \B) is a coherent superposition of pairs of 
particles with equal and opposite momentum (see Fig. [5]). 
Despite the strong similarity of Eq. (|36|) with the bound- 
ary states studied by Ghoshal and Zamolodchikov^ in 
the context of integrable scattering theories, however no- 
tice that integrability in scattering theories imposes some 
extra conditions on the kernel IC(k) which are not neces- 
sarily satisfied in the present case. 



III. HEISENBERG EVOLUTION AND 
TWO-POINT CORRELATION FUNCTIONS 

Let us now turn to the dynamics. In the lattice formu- 
lation, this can be studied by looking at the Heisenberg 
evolution^ of the fermionic operators Cfc(t). Setting 

c k (t) = u k (t)A k -v* k (t)(A°_ k y , (37) 

where A® are the eigenmodes of the initial Hamiltonian 
"H(ro), and using the fact that the time evolution of the 
operators Ak diagonalizing the final Hamiltonian is triv- 
ial, it is easy to obtain that 

c k (t) = Rt(T) u fc (t) R fe (r) R fe (r ) A° k , 



where 



U fc (t) 



-ie k t 







ie k t 



(38) 



(39) 



The relevant dynamical observables can then be evalu- 
ated on the time evolved state, Eq. by using the 
Heisenberg picture and Eq. (|37p . after writing them in 
terms of Jordan- Wigner fermions. We are now ready to 
provide some technical details on the computation of the 
two-point spin correlation functions. 



A. Order-parameter spin correlation functions 

The time-dependent correlation function of the order 
parameter is defined as: 



p^(t,r) = (B\al +r (t +t)o-Ut )\B} 



(40) 
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where to denotes the waiting time after the quench. Due 
to the translational invariance of the system, Eq. (|4T))) 
does not depend on m, and we can choose m = 1 with- 
out loss of generality. The operator a^ +r (to + t) cr^(t ) 
connects states with different c-fcrmionic parity; there- 
fore it cannot be simply evaluated using Jordan- Wigner 
fermions in the even Hamiltonian sector % + . This prob- 
lem can be circumvented^ by considering a four-spin cor- 
relation function on a chain of length L, with r < L/2: 

C*°(t, r; L) = (B\ *t_ r+1 (t + t) of (to) 

xal +1 (t a +t)al_ r+1 (t )\B) . 

By using the cluster property and taking the thermody- 
namic limits, one can show that this function reduces to 
the square of the correlation function p xx (t): 

{f*(t,r)) 2 = lim C xx (t,r;L); (42) 

the crucial advantage of this strategy is that the four- 
point correlator in Eq. (|41[) conserves the c-fermion par- 



ity, and can therefore be evaluated in the (antiperiodic) 
even fermionic sector—. 

Following Refs. [HlHoj], Eq. (j2TJ can be written as a 
Pfaffian. We first write it in terms of Jordan- Wigner 
fermions, by simply using a x = e 47r ^ fc <j " fc (cj + Cj); then, 
after defining Aj(t) = ct(i) + Cj(t) and Bj(t) = c\(t) - 
Cj(t), we get: 

C xx (t, r- L) = (B\ [BL +1 (t +*)••• B L _ r {t + t)} 
x [AL +2 (t Q +*)••■ A L _ r+1 (t + t)] 

x [Br (t ) ■ ■ ■ BL_ r (t )} x [A 2 (t ) ■ ■ ■ A±_ r+1 (t )} \B) , 

(43) 

in which we also used the equalities (1 — 2rij) = AjBj = 
—BjAj and {Aj,B t } = Vj, I (here {•,•} denotes the 
anticommutator) . An application of Wick's theorem then 
leads to the following expression, where the square of the 
correlator C xx is written in terms of the determinant of 
a (21/ - 4r) x (2L - 4r) matrix: 



(C xx (t,r;L)) 2 = 

(B Ul] (to + t)B [h] (to + t)) (B [n] (to + t)A [l2] (to + t)) (B [jl} (t„ + t)B [h] (to)) {B [n] (t Q + t)A [u] (to)) 

- (B [h] (to + t)A [n] (to + *)) (A [n] (to + t)A [h] (t + *)) (A []2] (t + t)B [h] (to)) (A [h] (t + t)A [h] (to)) 

- (B [h] (to + t)B [n] (to)) - (A [h] (to + t)B M (t )) (B M (t )B [h] (to)) (B w (t )A [h] (to)) 

- (B^to + t)A [M] (t )) - (A [h] (t + t)A [M] (t )) - (B M (to)A [u] (t )) (A [u] (to)A [u] (t )) 



' (Aj(ti)Ai(t 2 )) = ^e*U- l ^u k (h) + v k (t 1 ))(ut(1*)+vi(t 2 )) 

k 
k 

(Bj(ti)Ai(t 2 )) = i^e^- i )( Ufe (t 1 )-u fc (t 1 ))K(t 2 ) + ^(t 2 )) 

k 

(Bj(ti)Bi(t 2 )) =jJ2 e iHj ~ l) - («*(*2) ~ <(t 2 )) 



with < 



(44) 



(45) 



In Eq. (|44|) the two indexes j, I have subscripts which in- 
dicate their range: an index with label i runs from 4 + 1 
to L — r, with 2 runs from 4 + 2 to L — r + 1, with 3 
runs from 1 to ^ — r, with 4 runs from 2 to 4 — r + 1. 
Each of the 16 submatrices in (|4"4")l is therefore of dimen- 
sion (4 — r) x (4 — r), and can be explicitly evaluated 
by substituting the proper definitions of the Aj's and 
Bj's and then working in momentum space, as detailed 
in Eqs. (|45[) [fc-sums are taken for fermion antiperiodic 
boundary conditions, Eq. 0]. All the diagonal entries 
of the matrix in (l4"4")l are zero, since they do not enter the 
contractions in Wick's theorem 55 . 



B. Transverse-field spin correlation functions 

We now briefly sketch how to compute the two-point 
correlation functions of the transverse-field magnetiza- 
tion: 

p**(t,r) = (B\a* m+r (t +«)<(*<>) \ B ) ■ ( 46 ) 

As before, since we are using periodic boundary condi- 
tions, we can take m = 1. Contrary to the order parame- 
ter correlation function, the operator <J* l+r (t + 1) cr^(t ) 
conserves the c-fermion parity, therefore one can rewrite 
it using Jordan- Wigner fermions in the even Hamilto- 
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nian sector H + , and the computation becomes simple. 
Using a] = 2 c]c 3 - 1, one gets p zz (t, r) = (B\ (2c+ +1 (t + 

t)cr +1 (t +t)-l) x (2cJ(t )ci(t )-l) Switching to 
momentum space, and using Eq. (|37j) together with the 
fact that \B) is the vacuum for A k particles, the following 
expression can be easily obtained: 



^(t,r) = l-|E[K(*o)| 2 + k(*o + *)| i 



L 2 



£ [ e " ( '" fc) M*o + t) v* k (t ) Ul (t + t) u](t ) 

k,l 



+v k (t +t) vf(t )ui(t + t) u* k {t )) 
+ \Mto + t)\ 2 \vi(t a )\ 2 ^ 



C. Density of kinks 



(47) 



The density of kinks for the Ising model ([3]) is defined 
by the operator 



E 



(48) 



Its expectation value p^(t) at a certain time t on the 
boundary state \B) is strictly related to the nearest- 
neighbor order-parameter correlation function. For pe- 
riodic boundary conditions this can be simply written 
as: 



A*) = 



\-{B\aUtX+x{t)\B) l-p ra (0,l) 



(49) 

where p xx (Q, 1) is evaluated after a waiting time to = t. 
Proceeding in an analogous way as in Sec. IIII AIITTTBl 
this can be explicitly rewritten as: 



A*) 



-i (ujt (t)w; (i) - u£ (t)« fc (t)) sin fc] . (50) 



Notice that Eqs. (gSJ, 07} and (J3DJ) refer to a finite sys- 
tem of i spins; the thermodynamic limit can be formally 
attained by taking the continuum limit and replacing the 
sums with integrals: — * fdk/(2ir). 



IV. THERMAL RELAXATION OF NON-LOCAL 
OPERATORS: ORDER PARAMETER 
CORRELATIONS 

In this section we discuss in details the behavior of 
the correlation function of the order parameter Pq{t,r) 
for a quantum quench of the transverse field from To to 
r, focusing on the asymptotic relaxation for long times 



(t —> oo) and for long distances (r — > oo). The correlation 
function in Eq. (|40[) can be explicitly written as 



This correlator can be calculated following the prescrip- 
tions given in Sec. IIII A1 i.e., by numerically computing 
the Toeplitz determinant in Eq. (l44l) . 



A. Time-dependent correlations 

Let us start our analysis with the autocorrelation func- 
tions (r = 0). For the system at equilibrium, time trans- 
lation invariance implies that the autocorrelation func- 
tions do not depend on the waiting time to in Eq. (|5ip . 
Quite remarkably, we found that, as long as one con- 
siders pQ : (t,0) at large times, this is true also for the 
nonequilibrium case. Therefore, in this subsection we set 
for simplicity to = 0. The relevant features emerging for 
the case r = 0, to = have been already elucidated by 
some of us in Ref . (3~0 | . Here we summarize and discuss 
in more details the results anticipated in that paper. 

The envelope of the autocorrelation function 
Pq : (t,0) = (S|cr^(t)(T^ l (0)|K) relaxes exponentially 
to zero in time: 



'(t,Q) 



-t/T 



(52) 



This applies both to quenches (Fo ^ T) ending in the 
ferromagnetic phase (r < 1) or, as already pointed out 
by Calabrese and Cardy-ii, at criticality (T = 1). As it 
can be clearly seen from the left panel of Fig. [3J as soon 
as r =/= To, a sharp contrast with the zero-temperature 
behavior in equilibrium becomes evident (black curve). 
Indeed, at zero temperature and at equilibrium, the cor- 
relator PxLo(tj r ) asymptotically tends, both in time and 
in space, to the square of the order parameter (cr x )ii. 
In that case, the asymptotic value (a x ) 2 is not reached 
exponentially, but through an oscillatory power-law de- 
cay with an exponent a depending on the phase of the 
system (a = 1 in the ferromagnetic phase, a = 1/4 at 
criticality, and a = 1/2 in the paramagnetic phased). 

Drawing a parallel with the equilibrium physics, the 
exponential decay to zero observed for quantum quenches 
is reminiscent of the exponential decay of the auto- 
correlation function observed in equilibrium at finite tem- 
peratures T > 0. The rate Tq can be identified with the 
non-equilibrium analogue of the phase coherence time. 
In the equilibrium case, a semiclassical analysis due to 
Sachdev and Young^S shows that, for T -C A and in the 
ferromagnetic phase, one has: 



p x T x (t,r)^N 2 R T (t,r) 



(53) 



where Nq is the vacuum expectation value of the order 
parameter Nq — (a x ) and Rx(t,r) is a relaxation func- 
tion describing thermal excitations: 



Rr{t, r) = exp 



v k t 



(54) 
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Figure 3: (color online). Absolute value of the on-site time 
dependent correlation function pQ (t, 0) for quenches termi- 
nating in the ferromagnetic (F = 0.5 -left panel), and in the 
paramagnetic phase (F = 1.25 -right panel). Fhe various 
curves are obtained for different values of Fo, as indicated in 
the graphs. Fhe black curve corresponds to the equilibrium 
case (in the left panel, the horizontal dotted line indicates 
the asymptotic value {c x } 2 , while the dashed one denotes the 
power-law envelope ~ t -1 ). Colored curves stand for differ- 
ent quenches To ^ V; the straight dashed lines indicate the 
respective leading exponential decays. In the paramagnetic 
case, curves are rescaled by the zero-temperature 1 / yt pref- 
actor, analogous to the Kit, 0) term in Eq. (|58[) . 



Here Vk = dk£k is the velocity of the thermally ex- 
cited quasiparticles (kinks) and e _efc / T their Boltzmann 
weight. From Eq. (f54|) . one can readily extract the basic 
time- and length-scales of the system in the ferromagnetic 
region, the phase coherence time Tj, and the correlation 
length Indeed, R T (t,0) = e"*^ with: 

4=(/fe-^k|)"«^ e ^, (55) 

and Rx(0,r) = e~ r ^T with 

Z-Ut^T-M?"' (56) 

where c is the "speed of light" , given in terms of the band- 
width J and of the lattice constant a by c = 2 J a y/T. 
Notice that in the critical region (T ~ A) quasiparticles 
are not well denned and semiclassical arguments are not 
applicable. However, on the basis of an analysis of the 
continuum scaling limit one obtains 

t$ , & ~ ^ ■ (57) 

This behavior is confirmed by an exact analytic calcula- 
tion on the lattice^, giving Tji = -Ik. 

The parallel between equilibrium behavior of the au- 
tocorrelation function at finite T, and the behavior af- 
ter a quench is even more astounding when we consider 



quenches ending in the paramagnetic phase (r > 1), 
where the equilibrium analysis (in spite of the fact that 
quasiparticles are no longer kinks) predicts a very similar 
structure^: 

p^(t,r)~K(t,r)-R T (t,r), (58) 

with Rx(t,r) still given by Eq. ([Ml) , while K(r,t) is an 
oscillatory decaying function determined by the quantum 
fluctuations in the ground state. For the Ising chain, 
K(t,r) ~ Ko(A\/r 2 — Tt 2 /c), where K is the modified 
Bessel function, which oscillates and decays as 1/ for 
large t, while it decays exponentially as e _rA//c for large 
r. The data for the autocorrelator pq x {t, 0) as a function 
of time with final T > 1 show pronounced oscillations as 
in the equilibrium case, and the exponential drop clearly 
emerges only if pq x (t, 0) is multiplied by \ft, the leading 
decay of K(t, 0) (see the right panel of Fig. [3j) . 

Given a generic quench from Tq to T, one can numer- 
ically extract a relaxation rate Tq. As shown in Fig. 2 of 
Ref. [3l|, Tq decreases as the quench "strength" \Tq — T\ 
increases, implying that the more the system is driven 
out of equilibrium, the faster correlations decay in time. 
In the limiting case of equilibrium (r — > T ), Tq — > oo 
and the exponential relaxation turns into a power-law. 
Notice that, even though the decay rate of correlations 
depends on the strength of the quench, the value reached 

asymptotically is always zero, pQ X (t, 0) — 0, irrespec- 
tive of the values of Tq and T. 

We can understand these results as follows: the ef- 
fect of a quench from Tq to T consists of injecting in 
the system an extensive amount of energy. As seen in 
Sec IIII1 this leads to the generation of a finite density 
of quasi-particle excitations Ak , with a non-equilibrium 
occupation 

f k = {B\A\A k \B) . (59) 

The farther is Tq from T, and therefore the greater is the 
injected energy, the higher is the density of quasiparticles 
created. This is qualitatively very similar to the physics 
behind the decay of the autocorrelation function at equi- 
librium, if we consider the fact that the greater is the 
temperature T, the higher is the density of quasiparticles 
in the system. It is therefore quite natural to describe the 
dependence of the relaxation rate on the initial state and 
on the final Hamiltonian in terms of an effective tem- 
perature^. The most natural way to define the effective 
temperature T e g consists in comparing the energy of the 
initial state (B\H(T) \B) to that of a fictitious thermal 
state relative to the final Hamiltonian "H(T) 

(^(r)>T cff = ]Te£(n fe (T off ) + n_ k (T cS ) - l) , (60) 
fc>o 

where n/-(T e g) = (l + e e fe/ Toff ) _1 is the Fermi distribution 
function for the quasiparticles Ak diagonalizing H(T) [see 
Eq. (fT5|) ]. The effective temperature is therefore deter- 
mined by the implicit equation: 

(B\H(T)\B) = (H(r)) TQH , (61) 
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or, in a completely equivalent way, by: 

5>£/ fc = ^ e £ nfe(Tcff) _ (62) 

k>0 k><3 

For given To and T this equation always admits a single 
solution. 

Quite remarkably, Tq appears to be univocally deter- 
mined by T c g and by the energy gap A(r) in the final 
state. This means that two different quenches starting 
and ending at different Tq and T, but having the same 
effective temperature and the same A(r) will show the 
same Tq , irrespective of whether the two final T are both 
in the same region of the phase diagram or, for example, 
one in the paramagnetic and the other in the ferromag- 
netic region. Thus the fine details of the initial conditions 
are not important, while the only key physical parameter 
is the effective temperature of the initial state. 

As one can see from the main panel of Fig. |4l away 
from criticality the phase coherence time Tq (symbols) 
nicely follows, qualitatively and even quantitatively, the 
curve describing the equilibrium phase coherence time 
t t=t ft a ^ the effective temperature T e g [see Eq. (|55|) ] . 
This statement is true also at the critical point, T = 1, 
where, however, a finite-size scaling is mandatory in or- 
der to capture the T~g behavior predicted at equilibrium. 
Indeed, while the comparison of Tq with t^ =t ff system- 
atically improves as T c g is decreased, the long wavelength 
modes become increasingly important at low effective 
temperatures, thus making boundary effects more visi- 
ble. Numerical data are shown in the inset of Fig. 2] no- 
tably, after a finite-size scaling, the data not only follow 
the T~fi law but also appear to have the same prefactor 
8/ir analytically predicted in Ref. [5l| for a system at 
equilibrium. 

Since the parallel between quench dynamics and ther- 
mal behavior works better at low effective temperatures, 
one might be tempted to fix the effective temperature 
directly in the long-wavelength scaling limit. However, 
particular care should be taken in this case. Indeed, the 
energy of this state is easily computed using the matrix 
elements (|33[) and is given by 

h B = (B\U\B) = L / ^E(k)V 2 (k), (63) 
Jbz 27r 

where the explicit length L of the system here shows 
that we are in presence of a global quench. The inte- 
gral above is finite as long as the lattice spacing a is fi- 
nite, whereas it is logarithmic divergent as log(l/a) when 
a — > 0. The simplest way to see the appearance of this 
logarithmic divergence when a — > is to observe that, 
in the continuum limit, the mass couples to the operator 
i : ip+(x)ip-(x) : normal ordered with respect to the mass 
?7i of the field. If we start with too and change the mass 
m — > to = to + 5m, this operator acquires a vacuum 
expectation value, and the energy density per unit length 
changes as 

7<5to (B\tP + ( x )^( x )\B) = (5m) 2 | {p2 / A2)3/2 




Figure 4: (color online). Phase coherence time as a func- 
tion of the effective temperature, in the off-critical region for 
different values of F (main panel) and at criticality (inset). 
Symbols refer to Tq (here, as in Fig. [T] empty symbols stand 
for quenches starting with Tq < F; filled ones are for To > T), 
while continuous curves denote the equilibrium values at finite 
temperatures t!£, Eq. (|55p. For the system at criticality, we 
performed a finite-size scaling of Tq to study the asymptotic 
agreement with the equilibrium law t~j< = ~§T (bl ac k line), 
[data for V = 0.5, 0.75, 1, 1.25 acknowledged from Fig. 3 of 
Ref. HH]. 



This quantity, which is logarithmic divergent, matches 
with the limiting expression for a — > of the integral (|63l) , 
expanded at the same order in 5m. This implies that, 
keeping fixed the energy density per unit length h B / L of 
the boundary state while going to the continuum limit 
a — > 0, one is forced to consider only quenches where the 
differences in the masses are logarithmically close each 
other, to — to ~ l/log(l/a). With this adjustment of 
the masses, the physics of the quench process is invariant 
under a change of the lattice spacing. 

B. Space-dependent correlations 

A similar analysis can be performed for the space 
dependence of the order parameter correlation func- 
tions m, Pq (M = (B\a* l+r (to)a? n (t )\B), for suf- 
ficiently large times to after the quantum quench. As 
for the time autocorrelations, we find that this correlator 
displays an exponential decay to zero as a function of the 
distance r: 

^(0,r)^(r)e^. (64) 

This allows us to define the correlation length £q for a 
system after the quench, which is in general a function 
of the quenching parameters Tq and T. 

The behavior of Pq^O, r) as a function of the distance 
is explicitly shown in Fig. [5] For quenches in the ferro- 



11 





Figure 5: (color online). Space-dependent correlation func- 
tion PQ C (0,r) for quenches in the ferromagnetic phase (left 
panel, F = 0.5) and in the paramagnetic phase (right panel, 
r = 1.25). The various data sets are for different values of To, 
as explicitly indicated near each curve. On the right panel, 
dashed colored lines are exponential fits of the curves in the 
large-r limit, while the dotted-dashed black curve represents 
IC(r) at the corresponding T. In the insets we plot the cor- 
relation length £q as a function of To, as extracted from the 
exponential decay in the space of correlations. 



Figure 6: (color online) Distance r* at which the first max- 
imum in the oscillations for the spatial correlations crossing 
the critical point from the disordered phase appears, as a func- 
tion of F < 1 and for fixed Fo > 1. Such value is a measure 
of the period of oscillations. In the inset we show PQ E (r) for a 
paramagnetic-to-ferromagnetic quench. The oscillations are 
superimposed to an exponential decay behaving thermally, 
which is depicted with a straight red line. 



magnetic phase, one clearly identifies a pure exponential 
decay with r, as illustrated in the left panel of the fig- 
ure; therefore, the decay rate £q can be easily extracted 
by fitting the numerical data£2.. This is analogous to the 
equilibrium case, Eq. as already observed for the 

time-dependent autocorrelation functions. In the same 
way as for Tq, the correlation length increases as 
the quench strength decreases and eventually diverges for 
r = Tq (equilibrium, at T = 0), where at large distances 
Pt=o(0, r) tends to (<J X ) 2 and not to zero; the dependence 
of £q on Tq, for fixed T, is shown in the inset of Fig. [5J 

In the paramagnetic phase the situation is more sub- 
tle, as one can see from the right panel of Fig. [5] Here 
the behavior of p^f(0, r) at short distances is dictated by 
the prefactor K(r) which decays exponentially to a con- 
stant C 7^ 0, K(r) s» C + exp[— Ar/c] (a similar behavior 
is present in equilibrium, where K(0,r) ~ Kq(At/c), 
though with C = 0). For distances r ^> c/A, where 
K « C, the observed exponential decay univocally de- 
fines £q. As a function of T ff and A, this appears to 
qualitatively follow the law in Eq. (|56p , as in equilibrium. 

A similar peculiarity associated to the prefactor K(r) 
appears when the quench crosses the critical point from 
the paramagnetic to the ferromagnetic phase. Here the 
exponential decay of spatial correlators is superimposed 
to oscillations periodically changing the sign of the cor- 
relations, K(r) ss cos(r/r*). As shown in Fig.[6l the spa- 
tial period r* of these oscillations varies with the distance 
from criticality, and eventually appears to diverge at r c . 
We numerically verified that r* diverges for T — ¥ T c as 



r* ~ 1/ yT — 1, irrespective of the value of Fo (in anal- 
ogy to what analytically observed-^ for Tq — +oo). The 
presence of these oscillations has been first reported in 
Ref. 15] for a specific quench starting at Tq = +oo, and 
appears to be closely analogous to a similar phenomenon 
observed for linear quenches of the transverse field across 
the critical points. We will come back to this point in 
Sec. EH where we will show that these oscillations are an 
effect of the integrability of the model and that a suffi- 
ciently strong non-integrable perturbation of the system 
dynamics leads to their suppression. Despite these os- 
cillations, the leading exponential decay characteristic of 
thermal-like behavior is still present. Indeed, in anal- 
ogy with the phase coherence time, we find that the cor- 
relation length £q qualitatively follows the behavior of 
the corresponding equilibrium length ^ at the effective 
temperature T e g determined by Eq. (pjlj) . This can be ex- 
plicitly seen in Fig. where the values £g for quantum 
quenches are indicated by symbols, while the equilibrium 
data are plotted as continuous curves (see Sec. lIV Cl below 
for a discussion of the small discrepancies observed) . The 
dynamics at the critical point deserves a separate treat- 
ment, since the semiclassical analysis fails there. Again, 
after a finite-size scaling of numerical data, the data ap- 
pear to converge at low T e g to the behavior, as 
predicted for the equilibrium case at finite temperature, 
Eq. (|57|> . This is shown in the inset of Fig. [7] 
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Figure 7: (color online). Correlation length as a function of 
the effective temperature, in the off-critical region for different 
values of V (main panel) and at criticality (inset). Symbols 
refer to the quenched situation £gi continuous curves denote 
the equilibrium values at finite temperatures Eq. (|56|l . 
while dotted lines are obtained with a semiclassical analysis 
for the system out of equilibrium, Eq. ((65}. In the inset we 
show a finite-size scaling of data at criticality; the straight 
black curve indicates a behavior, and it is plotted as a 
guideline. 



Figure 8: (color online). Phase coherence times extracted 
by means of the semiclassical argument; data are shown for 
two values of V in the ferromagnetic phase. The continuous 
red curves indicate which is computed using the effective 
quasiparticle distribution nk, while the dotted blue ones are 
for Tq and are computed according to the occupation factor 
fk- We also show the phase coherence times Tq evaluated 
numerically from the non-equilibrium situations (symbols - 
same data of Fig. |4}. In the inset we plot rik and fk for a 
particular value of (Fo = 0.75, T — 0.5), corresponding to the 
point in the main panel that is indicated with an arrow. 



C. Semiclassical analysis out of equilibrium 

A better agreement between numerics and the semi- 
classical analysis can be obtained by modifying the anal- 
ysis of Ref. [40] for the quench case, by substituting the 
Boltzmann weight e -Cfc / T with the occupation factor fk 
for the eigenmodes after the quench. In analogy with 
Eq. (|56|) . we obtain 



<p _ 



dk. 
— Jk 

7T 



(65) 



Results of the computation using this formula are shown 
as dashed curves in Fig. and reveal a marked improve- 
ment in the quantitative agreement between numerics 
and theory. A similar analysis applied to the autocor- 
relation time leads to a generalization of Eq. ([55)1 into: 



(— f nfe(T)|vfe|(ifc) _1 . This gives a result that essentially 
coincides with the formula for given in Eq. (1551) . The 
two phase coherence times Tq and in the low temper- 
ature regime are very close, apart from some constant 
prefactor, in spite of the fact that the two quasiparticle 
distribution functions fk and rik can be rather different, 
as explicitly shown in the inset of Fig. |8j for a specific 
value of (r ,r). 

Finally, we mention that a double check of the con- 
sistency of our analysis comes from the comparison be- 
tween the effective temperature evaluated from Eq. (|61[) 



and that obtained using the equality Tq = fj,. In this 



last case, in the limit T c g <C A one finds 



T, 



A 



cir 



in[A/(r - r ) 2 ] 



(67) 



dk 

— Jk\Vk 
7T 



(66) 



Data displaying Tq for two values of T are shown in Fig. [8] 
(blue curves). As one can see, the small discrepancies be- 
tween Tq and r^ =T(sff at low temperatures that appeared 
in Fig. 3] are now reconciled within numerical accuracy. 

We point out that one could also evaluate the phase 
coherence time by using the effective thermal quasipar- 
ticle distribution function nk(T c g) in Eq. (|66[) . instead 
of the non-equilibrium distribution function f%,= 



in agreement with the cusp that emerges for T e fj as one 
approaches the limit r — > T (see the inset of Fig. (2) in 
Ref. M)- 



V. NON-THERMAL BEHAVIOR OF LOCAL 
OPERATORS 

In this section we focus on a different class of oper- 
ators, belonging to the local sector with respect to the 
fermionic quasiparticles. We first consider the two-point 
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time-dependent correlation functions of the transverse- 
field magnetization tr z , and then discuss the behavior of 
the nearest-neighbor correlations of the order parameter 
a x , which is equivalent to the density of kinks. 



A. The transverse-field spin correlation functions 

The two-point time-dependent correlation functions of 
the transverse-held magnetization a z are defined as: 



p* z (t,r) = (B\** m+r (t +t)a z m (t )\B) 



(68) 



in a way analogous to the order-parameter correlation 
functions, Eq. (pfU|). These can be evaluated as sketched 
in Sec. lIIIBl As for the case of Pq^O, r), we will consider 
a waiting time to long enough that results are indepen- 
dent of it. 

The correlation functions of the transverse-field opera- 
tor a z for the system at equilibrium (To = T) have been 
computed in Ref . [l| , where it has been shown that they 
relax to the square of the transverse magnetization (er z ). 
This is given, in the thermodynamic limit, by 



(a z ) T = — 

IT 



cos(2Afc) tanh 



(-) 

\2TJ 



dk 



(69) 



in such a way as to have 



with A fc = ^arctan( sil } k r ) 
< Afe < 7r. In contrast to the correlations of the order 
parameter, at equilibrium p%?(t, r) exhibits a character- 
istic power-law decay in time, both at zero and at finite 
temperature: 

\pl?(t,r)-(**)l\~t-«T , (70) 

which makes it impossible to define a time-scale analo- 
gous to Tq . The power- law exponent ax depends on the 
phase of the system and on the temperature. We cal- 
culated «t by numerically fitting the curves of p^ z (t,r) 
(evaluated from exact analytic formulas- 1 ) for long times 
t. At zero temperature we found: 



a 




for r < 1 
for r = 1 
for r > 1 



(71) 



On the other hand, at finite temperatures ar = 1, irre- 
spective of the system phase. More precisely, we found 
that the ay = 1 behavior always holds for long times, 
while at short times and for T > 1 correlations decay 
with an exponent 2, up to a given transient time t* which 
decreases with the temperature. 

We now consider the correlations of the transverse field 
p z q{t, r) in the case of a quench and show that, contrary 
to the behavior of the order parameter correlation func- 
tions (foTj). they do not exhibit thermalization. For this 
purpose, we first concentrate on the asymptotic values 
that are reached at long times or distances, and then 
focus on the exponent of the power-law temporal decay. 
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Figure 9: (color online). Asymptotic value reached by the 
transverse- field correlation functions, which is the square of 
the magnetization (u z ). The magnetization for a quenched 
dynamics (symbols) corresponding to a given effective tem- 
perature, as defined in Eq. (|61|) . is compared to the corre- 
sponding value at equilibrium, at the same temperature (con- 
tinuous lines). Different colors and symbols stand for the 
various values of F, as declared in the legend: F = 0.5 (black 
circles), 0.75 (red squares), 1 (green diamonds), 1.25 (blue tri- 
angles). Inset: {<j z )q reached asymptotically after a quench 
from a given r , to T = 0. Notice that the corresponding 
equilibrium value {<J Z ) T for T = is rigorously zero, at any 
temperature. 



As in the equilibrium case, the two-point transverse 
spin correlation functions for the system out of equilib- 
rium relax asymptotically, both in time and in space, 
to the square of the expectation value of the transverse 
spontaneous magnetization, provided the waiting time to 
after the quench is sufficiently large: 



p zz (t,r) 



Contrary to the magnetization along the coupling direc- 
tion (a x ) T , which is rigorously zero at any finite temper- 
ature, the transverse magnetization {<J Z ) T is finite even 
at T > 0, see Eq. It would then be natural to 

carry out an analysis on (<J Z ) Q , analogous to the one per- 
formed for the relaxation dynamics of the order param- 
eter correlations, to check if it has anything to do with 
the equilibrium transverse magnetization (c z )t c!! - 

We have found, however, a completely different behav- 
ior. A signature of this fact is evident from the inset of 
Fig. [9l where we plotted the value of the transverse mag- 
netization (<j z ) q which is asymptotically reached after a 
quench towards T = 0. As one can clearly see, contrary to 
the finite-temperature equilibrium value (c z ) T which is 
rigorously zero for T = at any finite temperature T [see 
Eq. (|69l) ]. the value after the quench is always non-zero, 
as long as To 7^ 0. To be more quantitative, after defin- 
ing an effective temperature for the out-of-equilibrium 
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system according to Eq. (|6Tj) . one discovers that 

Ms * (Or*, • 

This is explicitly shown in Fig. El where we plot (<r z ) Q as 
a function of T e g (filled symbols), together with the cor- 
responding values (<t z ) t at the same temperature (con- 
tinuous lines). Contrary to what is observed in Figs.0Kl 
here we do not find a thermal behavior, except for the 
very specific case where the system is quenched towards 
the critical point T c (green diamonds and line in Fig. [9J . 
This can be understood in terms of the behavior of the 
quasi-primary operators of the Conformal Field Theory 
of the critical point^. 

Let us now briefly discuss the finite-time transient 
behavior of Pq (t, 0) for the system out of equilibrium. 
In a way analogous to the finite-temperature case, the 
transverse- field correlations p Z Q Z (t,Q) relax in time as a 
power-law: 



r 



We found the following power-law exponents olq, depend- 
ing on the system phases before and after the quench: 



a Q 


r < i 


r = i 


To > 1 


r < 1 


i 


i 


2 


r = i 


2 


3/2 


2 


r > i 


i 


1 


2 



Since for the system at equilibrium and at any finite tem- 
perature we have ax = 1, this means that, in general, 
even the finite-time power-law behavior is a non-thermal 
one. The jump in the values of ao in Eq. (fTTj) and of a Q 
in the table, depending on the different system phases 
before and after the quench, can be seen as due to lattice 
effects. 



B. The density of kinks 

As a further confirmation of the non-thermal behav- 
ior of local operators, we concentrate on the density of 
kinks p 1 ^ . In particular we compare the asymptotic value 
reached after a quench (p^)q with the thermal behavior 



at equilibrium. 
byi2^: 



In the thermodynamic limit this is given 



(P 



dk 
2t7 



1 - 



1 + T cos k 



vi + r 2 + 2r cos k 



(l-2n k (T)) 



On the other hand, the out-of-equilibrium behavior of 
(p^) q = limt->oo p^ (t) after a quench is retrieved from 
Eq. (|50[) and by taking the asymptotic long-time limit for 
L — > co. 

The comparison between the values of {p^)q after a 
quench and {p^) T „ at equilibrium, where T e g for the 
corresponding out-of-equilibrium system is obtained from 
Eq. (|6ip. is presented in Fig.[TU] As one can see, outside 
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Figure 10: (color online). Asymptotic value of the density of 
kinks as a function of the effective temperature T e ft. Different 
colors stand for various values of F, as shown in the inset. 
Symbols denote the density of kinks {p^)q after a quench; to 
each value of To corresponds a different initial state, therefore 
a different effective temperature T e g . Straight lines denote the 
finite-temperature equilibrium values (p J ^)r=T off - 



criticality the two quantities are evidently not related 
and behave in different ways. On the other hand, the 
results for quenches ending at the critical point r c = 1 
(diamonds) perfectly follow the dotted-dashed equilib- 
rium line (green data) , so that the density of kinks after 
a quench to the critical point is univocally determined 
by the effective temperature T c ff. Remarkably such be- 
havior is not found for a non-critical dynamics, where 
fine details of the initial condition seem to be important. 
As explained above, we frame the effective thermaliza- 
tion only at criticality in terms of the behavior of the 
quasi-primary operators in conformal field theory^. 

Incidentally, the thermal behavior of {p^)q at critical- 
ity can be recovered analytically. Indeed, in that case the 
density of kinks at equilibrium with an effective temper- 
ature T e ff is given by substituting r = T c = 1 in Eq. (|72[) : 



' r cff I 



dk 


1 


2t7 





)e r k ° [l-2n fe (T eff )] 



(73) 



where e£ c = 2 a/2 — 2 cos k is the energy of the Ak quasi- 
particle at r c = 1, and the second equality follows from 
Eq. dHJ) with n k (T) = n_ fc (T). The energy {U(T)) Tea 
of the system after the quench is obtained directly from 
Eq. (|oT)]) . by evaluating (B\A\.Ak\B), and is given by 

r * dfc2(l + r )(l +cosfc) 



Mr)>T Bff 



On the other hand, the value of (p ) Q for a quench at 



(74) 
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criticality is extracted from Eq. (ISTfl) for f ^ oo: 

v , r*/i p l + rpcosfc o (r -l)sin 2 fc ) 

which, after simple algebra, can be shown to reduce to 
Eq. (|73j) . hence obeying the rigorous equality 

</0qL= ■ (75) 

The very similar behavior of (cr z ) and (p^) after a 
quench follows from the structure of the Ising Hamilto- 
nian in Eq. ([3]) , which is a sum of the kink-operator term 
plus a transverse magnetization term, and the fact that 
the total energy is a conserved quantity in the quench 
process. 

C. Description in terms of a generalized Gibbs 
ensemble 

The non-thermal behavior of the long-time asymptotic 
values for the transverse magnetization {o z ) Q and for 
the density of kinks (/A) Q outside criticality can be rec- 
onciled in terms of a generalized Gibbs ensemble pg 1 ®, 
which takes into account all the non-trivial constants of 
motion that typically prevent integrable systems from 
having an ergodic behavior: 

PQ = ^exp ( - V/3 fc e£n fc ) , (76) 

where Zg — Tr[e _ ^ ^ fc£ fc" fc ] is the partition function, 
and nk = A^Ak is the number operator of the Bogoli- 
ubov quasiparticle with energy e£. In Eq. (j76[) the op- 
erators {nk} assume the role of a full set of constants 
of motion, thus keeping track of the details of the initial 
condition. The fictitious inverse temperature fik is taken 
to depend on mode k of the system, and is a Lagrange 
multiplier that can be calculated from an equation for the 
expectation value of the corresponding integral of motion 
nk for the evolving system after the quench H(T), on the 
boundary state \B): 

(n k )g = (1 + eM)- 1 = (B\ n k \B) = f k ■ (77) 

After expressing the number operator n k in terms of the 
creation and annihilation operators for the initial Hamil- 
tonian H(To), one arrives at the following expression for 
the r.h.s. of Eq. ([77]): 

fk = (0(T O )| n k |0(T O )> = \ - ala }l ± bl , (78) 

L le k e k 

which, together with Eq. (|77p. provides an explicit for- 
mula for the /3fc 's^- 
Using such obtained fif., one finds that the asymp- 
totic values of both the transverse magnetization and the 



density of kinks after a quench are identical to the cor- 
responding expectation values in the same generalized 
Gibbs ensemble, irrespective of the initial and the final 
transverse fields r ,r: 

MsHOs (f^) Q = Wg. (79) 

This shows that the Gibbs distribution ([761 . with the 
same fictitious temperatures f3k defined in Eq. (j77[) . is 
able to catch the non-thermal behavior of both a z and 
Af after the relaxation outside criticality. As pointed out 
before, the two quantities necessarily have to behave in 
the same way, due to the energy conservation after the 
quench and to the particular structure of the Hamilto- 
nian ([3]). 

VI. SENSITIVITY TO THE BREAKING OF 
INTEGRABILITY 

In this section we briefly address the problem of how a 
non- integrable perturbation of the Ising model Q can af- 
fect the relaxation to the steady state. We do not intend 
to give here a fully comprehensive study of the effects 
of integrability breaking: more modestly, the purpose of 
this section is to provide evidence that some features of 
the asymptotic equal-time correlators, in particular the 
spatial oscillations after a quench, disappear as soon as 
the non-integrable perturbation is increased. This sheds 
light on the scenario we proposed throughout the pa- 
per: strong effects on the local observables (which do not 
thermalize by themselves, in integrable systems) are ex- 
pected, while the qualitative behavior of non-local ones 
(which, on the contrary, are already "thermal") should 
not change. 

One way to break integrability of the Ising model is to 
add a next-to-nearest-neighbor coupling of strength Ji : 

n 2 (T) = -jjr \a*a* +1 + IV*] - J 2 ]T of of +a . (80) 

For such non-integrable system, the standard analytic 
tools based on the Jordan- Wigner transformation fail, 
therefore one generally has to resort to a fully numeri- 
cal treatment. Here we employ an exact diagonalization 
study of model ([50)1 : unfortunately, this severely limits 
the actual sizes of the systems under consideration to 
- O(10) spins. 

We discuss the effective "generic" thermalization in- 
duced by integrability breaking, by analyzing the changes 
in the space asymptotics of the order parameter correla- 
tors, after a quench in the transverse-field strength T: as 
discussed in Sec. IIVB1 for quenches crossing the Ising 
quantum critical point from the paramagnetic to the fer- 
romagnetic side, non-thermal sinusoidal oscillations in 
real space naturally emerg o 15 ' 52 . Below we numerically 
show that a non-integrable perturbation leads to a sup- 
pression of them, thus driving the system towards a com- 
pletely thermal-like behavior, where only the exponential 
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Figure 11: (color online). Fourier transform of the equal-time 
order-parameter correlation function for the Ising model in 
presence of a non-integrable perturbation J2. Symbols con- 
nected by lines denote n^f(k) after a quench from To = 5 
to r = 0.4, according to the diagonal ensemble prediction. 
Dashed/dotted lines stand for the corresponding expectation 
values in the canonical ensemble n" ff (k) at the corresponding 
effective temperature T e g. Inset: absolute differences Sn xx (k) 
between the diagonal and the canonical ensemble predictions. 
Data are for a chain with L = 12 sites. 



decay in space of correlations is present. Instead of con- 
sidering the correlations ^^(O^) in real space, here we 
work in the momentum space, where the effects of the 
oscillations clearly emerge also for very small sizes. We 
define 



1 L 



(81) 



r=l 



In presence of a pure exponential decay in space of 
correlations p X Q ' , a typical Lorentzian curve for n xx (k) 
emerges. The superposition of a simple sinusoidal modu- 
lation induces the formation of two humps in fc-space, as 
it can be seen from the black continuous curve (circles) 
in Fig.HU for J 2 = 0. 

In Fig. [TT] we numerically compare the expectation 
value riT X ff (k) in the canonical ensemble at an effective 
temperature T e ff as given by Eq. (|61l) . with the asymp- 
totic value after the quench rvq(k), calculated from the 
diagonal ensemble^: 



n x Q x (k) = lim (4, t \n xx (k)\4>t) = 



Ei 



(<Pi\n**(k)\<pi) , 
(82) 

where \(j>t) is the state at time t after the quench \<j>t) = 
e -i« 2 (r)t |^ are .j- ne eigenstates of the Hamiltonian 
after the quench 'M.2(T), and c,; = (B\ipi). Humps are 
present in the out-of-equilibrium system for very small 
values of J2, while they gradually disappear when in- 
creasing the perturbation; on the contrary, they are com- 



pletely absent at equilibrium at finite temperature, thus 
clearly suggesting that they are typical non-thermal fea- 
tures. A clear evidence of the thermalization induced 
by non-integrability is given by the difference 8n xx (k) 



[k)\ as a function of J2. As shown in the 
(k) is decreasing with J2, until this 
becomes large and the system turns out to be again in- 
tegrable. 



inset of Fig. [TTJ Sn 



VII. CONCLUSIONS 

In summary, we have analyzed the temporal relaxation 
of some observables after a sudden quantum quench in 
a completely integrable one-dimensional spin model ex- 
hibiting a quantum phase transition. In particular, we 
focused on the spin- 1/2 quantum Ising chain, where we 
quenched the transverse magnetic field and studied the 
space and time asymptotics of two-point spin correlation 
functions corresponding to the spin operators along the 
longitudinal and the transverse direction, and the density 
of kinks. 

We found that, despite the complete integr ability of 
the Ising model, some observables exhibit thermal be- 
havior, while others do not. The thermal behavior is 
typical for operators which are non-local with respect 
to the fermion quasiparticles that diagonalize the model. 
The spin operator along the coupling direction, which 
is the order parameter, is of such type. In this paper 
we have explicitly demonstrated that correlation func- 
tions of the order parameter drop exponentially to zero 
and thermalize to an effective temperature which is ruled 
by the energy of the initial state after the quench. On 
the other hand, as examples of operators which are lo- 
cal in the quasiparticles, we considered the spin opera- 
tor along the transverse field direction, and the two-point 
nearest- neighbor correlator of the spin operator along the 
coupling direction. We showed that such two-point cor- 
relations present non-thermal power-law decay to some 
residual value, which is not related to its thermal counter- 
part at the associated effective temperature, except in the 
very specific case of a quench towards the critical point. 

It is worth stressing that after breaking the integra- 
bility of the model by adding an extra operator to the 
Hamiltonian [a next-to-nearest neighbor coupling, see 
Eq. (|80[)]. spurious non-thermal effects of non-local oper- 
ators disappear, while the leading thermal behavior does 
not qualitatively change. In particular, we have provided 
numerical evidence for the suppression of non-thermal si- 
nusoidal oscillations in real space of the order parameter 
correlators, after a quench from the paramagnetic to the 
ferromagnetic phase. 

We point out that the relaxation dynamics of a many 
body system following a quantum quench has been con- 
sidered also in other similar contexts. In Ref. [HH the 
behavior of the staggered transverse magnetic moment 
in a XXZ spin chain initialized in the perfect antiferro- 
magnetic state was analyzed: the staggered moment al- 
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ways decays exponentially, similar to our findings on the 
order parameter correlations in the Ising model. In some 
cases, depending on the value of the anisotropy param- 
eter, oscillations in the dynamics of the magnetization, 
which resemble the ones we reported for quenches from 
the paramagnet to the Ising ferromagne t 15 ' 52 , were also 
found. The time evolution of an initial state equivalent to 
the Neel state has been also studied for the Bose-Hubbard 
model^. In this case, oscillations of local observables 
such as the local density or nearest-neighbor correlators, 
seem to decay algebraically for all values of the interac- 
tion, while no crossover to a non oscillatory regime has 
been found. The differences between the two cases have 
been ascribed to the role played by the crossing of equi- 
librium critical points. 

In light of all these results, it would be clearly inter- 
esting to investigate other quantum integrable models, 
in order to see whether thermalization behavior occurs 



in correlation functions of operators that couple to in- 
finitely many quasi-particle states of the theory, while 
it is absent in correlation functions of operators with a 
finite number of matrix elements. 
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